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Abstract 

We report the result of our evaluation of the Feynman diagrams appearing 
in the determination of the four-loop renormalization group functions in the 
two-dimensional lattice O(n) cr-model by Caracciolo and Pelissetto. In the 
list of the integrals used for the determination of those functions, we hnd 
that three entries were not correctly evaluated. We give the values for them 
corrected by us including those for all other integrals which we computed with 
higher precision. These results are then applied to revise the determination of 
the second analytic correction to correlation length £ and spin susceptibility 
X by Caracciolo et al. as well as our determination of the mass gap by means 
of a hnite volume technique where we explicitly made use of the four-loop 
/3-function. In both cases we hnd sizeable changes in predictions. In the 
meantime there appeared a paper by Alles et al. where they revised one 
hnite integral in the list of our corrected integrals. After having taken the 
new revised value into consideration, we found that there are no noticeable 
changes in the perturbative predictions of the present paper including the 
hnal conclusions. 



1 Introduction 



The two-dimensional non-linear O(n) a- model is interesting in theoretical physics 
in many respects. In condensed matter physics, the model can be applied to the 
study of ferromagnetic systems. In elementary particle physics, on the other hand, 
it serves as a good toy model for the investigation of the strong interactions since it 
shares with QCD many common properties. 

The most interesting property of the a-model from the view of particle physicists 
is that of asymptotic freedom which is predicted by perturbation theory. Besides 
this, the model shows many additional interesting properties. One believes a non- 
perturbative generation of a mass gap via dimensional transmutation which controls 
the exponential decay of correlation functions at large distances. From the existence 
of an infinite set of conserved non-local quantum charges [[I], @, || |4j] , one can show 
that there is no particle production in the theory and also derive "factorization 
equations" directly. These properties lead finally to the exact determination of the 
S'-matrix up to a CDD-ambiguity [|, ^ , assuming that the spectrum of the theory 
is known. By making use of the exact S'-matrix and the thermodynamic Bethe 
ansatz, Hasenfratz and Niedermayer determined the ratio of the mass gap to the 
A-parameter m/A of the theory exactly 0, f|. If one formulates the a-model on the 
lattice, it has a further advantage: due to low dimensionality it is simple to simulate 
numerically and can therefore be used in many cases for testing new ideas in lattice 
theories. 

The property of asymptotic freedom is based on the validity of perturbation 
theory. One usually formulates perturbation theory on a non-symmetric trivial 
vacuum. The Mermin- Wagner theorem 0, however, tells us that in O(n) cx-model 
the continuous 0(n) symmetry can not be spontaneously broken in two dimensions. 
Therefore, the validity of the perturbation theory in infinite volume is in fact not 
guaranteed. This observation, together with the lack of clear asymptotic scaling in 
simulations, gives some authors 1TDL O, [E| a motive for an argument that the non- 
linear 0(n) a-model has no mass gap and for any n > 2 it undergoes a Kosterlitz- 
Thouless phase transition O, [Tj| at a finite coupling f c ^ 0. 



Monte Carlo simulations allow in principle to explore the behavior of a model in 
all areas of the bare coupling. It is of general interest to connect Monte Carlo results 
with those of the perturbative renormalization group (RG) which should dictate the 
scaling laws of the theory near the fixed point fo = 0. For a precise comparison 
with the Monte Carlo data, the computation of the RG functions at high orders 
in perturbation theory is desirable. Until now they are known up to four loops 
where the long and difficult four-loop computation was recently done by Caracciolo 



and Pelissetto [JL5|]. These determinations of the four-loop j3- and 7-functions then 
allow us to compute the correlation length £ and spin susceptibility \ up to second 
analytic correction. 

We applied an efficient coordinate space method for the evaluation of lattice 



Feynman diagrams described in detail in our paper Jl6| to compute the integrals 
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listed in the paper by Caracciolo et al. more precisely. In most integrals, our 
results showed good agreement with those by Caracciolo et al. We, however, found 
that three integrals were not correctly evaluated. Since these errors turn out to 
cause sizable changes in the perturbative prediction to £ and x as wen as m our 
determination of the mass gap by a finite volume technique jl7 l, we would like 
to report our results for the evaluation of the lattice integrals together with the 
corrections to the physical prediction to £ and \ as well as to our determination of 
the mass gap. Furthermore, the four-loop RG functions evaluated by Caracciolo et 
al. were applied in many papers [18|, [19|, |2(| [2l|, |12| and it is necessary to revise 
them in accordance with our corrected values for the RG functions. 

We would like to stress that we checked only a minor part of the four-loop RG 
functions by Caracciolo et al., where we just reevaluated the integrals listed in that 
paper. We still do not know whether the other analytic formulas are correct. For 
the complete confirmation of the four-loop RG functions, it would take a while until 
we finish our own computation of them which is in progress ||22|| . 

The paper is organized as follows. In the next section we define the cx-model 
and give in this model the hitherto known results for the perturbative RG functions 
including the four-loop determination by Caracciolo et al. and its corrected version 
by us. Section [3] is then devoted to the reinvestigation of the physical quantities 
by making use of our corrected values for the RG functions. Finally, the explicit 
evaluations of lattice Feynman diagrams are presented in the Appendix. 



2 Computation of the RG functions 



We consider the non-linear O(n) cr-model in two dimensions defined by the action 



S 



1 

2/o 



d x 



df,q{x) ■ d^q(x) - 2hq n (x) 



(2.1) 



where qi(x) with i = l,---,n is an n-component real scalar field of unit length 
q(x) 2 = 1 and fo denotes the bare coupling constant, h is here a magnetic field which 
regularizes the infrared divergences and explicitly breaks the 0(n) symmetry, lining 
up the field along the component q n . At small coupling, one can then consider small 
fluctuations around the direction of the magnetic field. In this way, the perturbative 
expansions in f Q can be obtained by setting 



q(x) = \iri(x), - ■ ■ ,ir n „i(x), 



ii—i 

a=l 



1/2 



(2.2) 



which solves the constraint q(x) 2 = 1. The regulator h can be removed at the 



final stage of the computation. This is guarenteed by the conjecture of Elitzur [23 



proved later by David order by order in perturbation theory ||24| , which states that 
in 0(n) cr-model, 0(n) invariant correlation functions have a finite infrared limit in 
two dimensions. 
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With this prescription, one can now start to compute the physically interesting 
quantities according to perturbation theory Higher loop computations are however 
accompanied by ultraviolet divergences, which requires a renormalization. This can 
be achieved by introducing a regularization and redefining bare fields ir and coupling 
/o as renormalized ones, 7Tr and /. In the dimensional regularization with dimension 
d, one sets n R = Z^ x I 2 tx and / = Z-f 1 f fi d ~ 2 , where the Z and Z 1; taken as a series 
in /, are adjusted so that the renormalized theory satisfies a chosen set of two 
independent conditions. Choosing \x as the scale of the renormalized theory, if one 
changes /i, one has to change the renormalized fields and coupling accordingly, which 
is represented by the RG equation 



4 +/3(/) ! + ^ t(/)+ 



\iU) + jP(f) 



h 



d_ 

dhi 



}G§W,/^) = 0, (2.3) 



where G$ is the renormalized n-point correlation function for the 7Tr field and /?(/) 
and 7(/) are (3- and 7-functions respectively: 



(3(f) = ^f = -b 1 f 2 -b 2 f 3 -b 3 f-hf + 0(f) 

7 (/) = ^\nZ = Cl f + c 2 f 2 + c 3 f + c 4 f + 0(f) 
d/i 



(2.4) 
(2.5) 



The coefficients bi, b 2 and C\ are universal in the sense that they do not depend on 
the renormalization scheme chosen and explicitly given by 



h 
b 2 



n 



2tt ' 
n-2 

2?r 



(2.6) 
(2.7) 



All other higher terms are, on the other hand, scheme dependent. In the minimal 



subtraction (MS) scheme, they were computed up to four loops |25], |2f| 

In 2 -A 



h 
c 2 
c 3 

C4 



4 (2tt) 3 
n-2 



12 



[n 



22n + 34) + -C(3)(n-3) 



32tt 3 
1 



1927T 



(n-l)(n-2), 

l)(n-2)[4(5-n) + 3(3-n)C(3)]. 



{11 



(2.9) 

(2.10) 
(2.11) 
(2.12) 

(2.13) 



with C(3) ~ 1.2020569. 
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If one formulates the cr-model on the lattice, there are many different ways to 
define actions to regularize the theory. Here, we would like to consider the standard 
nearest-neighbor action 



1 



Jo 



(2.14) 



The RG equation on the lattice which is formally the same as Eq. ( [2.3|) , but written 
in terms of the bare j3- and 7-functions, is given by 



Jt(/o)+|/3(/o) 
1 Jo 



h ^} G L\pJoA/a;h) = 0, (2.15) 



where G L is the lattice n-point correlation function for the ir field and the lattice 
(3- and 7-functions are given by 



7(/o) 



-a—h = -&i/o 2 - hfo - hfo - Wo + 0(JS), 
In Z = cjo + hfl + a 3 / 3 + cJo + Wo)- 



(2.16) 
(2.17) 



The coefficients in this regularization are also known to four loops. Up to three 
loops, they are given by |29|, |30|, |31 



61 

k 

h 

Cl 
C3 



&2, 

n-2 

Cl, 

n — 1 

8tt 
n — 1 



- + -7T 2 - 47T 2 G^ (n-2) + 1 + -IT - — 7T 2 

2 8 V V ; 2 24 



(27T) 



- - -7T 2 + 4tt 2 G 1 ) (n - 2) + — 7T 2 
2 8 V V ; 24 



(2.18) 
(2.19) 

(2.20) 

(2.21) 

(2.22) 

(2.23) 



with G\ ~ 0.0461636. Recently, Caracciolo and Pelissetto |15j computed the quite 
demanding four- loop coefficients 64 and C4: 



n - 2 f 2n - 7 n 2 - 4n + 5 



2tt 'V 96 32tt ' 16tt 2 

n — 2 . . n — 2 , , n — 6 
(n - 2 + 37r)Gi (3 - 7r)i? + —7— J 



+ 



77 



5n - 9 73n - 164 



16vr 3 



967T 2 



C(3) 



7T 



+ 



2tt 

+ i(n - 3)(n - 2){K - 2V 3 ) 

n-2 



24 



2Vi + (n - 2)V 2 + 41/ 4 - 2U 5 - 2V 6 - 16W 2 



(2.24) 



c 4 



n — 1 

2tt 
n-2 

~192tt 3 
n — 6 



5n - 21 (n - 2) 2 



3(n 



192 

(6n + 37)C(3) 
2 



32tt 
n-2 r 



47T 



327T 2 

2fn - 2) 



2) // 



167T 3 



5tt 



n 



2tt 



'■(3-n)R 



J 



n 



+ 



48 
n-2 



-Li 



(n - 3)(n - 2)(K - 2V 3 ] 



3 16 

[n - 3)^i + 2(n - 2)V 2 - AV A + 2V 5 + 2V 6 + 16W 2 



(2.25) 



where the integrals Gx, R, J, K, L±, V\, ■ ■ ■ , V§, W\ and W% are defined in Appendix ^ 
and their numerical values evaluated by Caracciolo et al. are listed in the middle 
column of Table |l[ 



Integrals 


Values I 


Values II 


G x 


0.0461636 


0.04616362923(1) 


R 


0.0148430 


0.014842965985741(1) 


J 


0.1366198 


0.136619772367581(1) 


K 


0.095887 


0.09588764425(1) 


Lx 


0.0029334 


0.002933289310(1) 


Vi 


0.016961 


0.016961078576(1) 


v 2 


-0.00114 


-0.0011381683(1) 




0.07243 


0.07243631946(1) 


v A 


-0.0013125 


-0.002624929(1) 


v 5 


0.010063 


0.0100630371(1) 


V 6 


0.017507 


0.017506946(1) 




-0.0296860 


-0.007421482994(1) 


w 2 


0.00221 


0.0006923019(1) 



Table 1: Numerical evaluations of the lattice integrals appearing in the 
four-loop (3- and 7-functions. The numbers in the middle column (Values I) 
are the evaluations by Caracciolo et al., while the numbers in the right 
column (Values II) are those by the author. 

By applying the coordinate space method for the evaluation of lattice Feynman 
diagrams described in detail in our paper [ IB] , we reevaluated all integrals in Table |l] 
with higher precision. We show our results in the right column of the same table. 
As we see in Table [I], we reproduce most values by Caracciolo et al. with a minor 
precision mistake for the constant L\. However, for three integrals^ V4, Wx and W 2 , 
our numbers are different from those by Caracciolo et al., where the difference in V4 
can be accounted for by a factor of 2. 

1 Recently, Alles et a l. |35| revised the value for W<i. We improved its numerical precision to get 
W 2 = 0.0021225468(1) " 
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We remark that Caracciolo et al. checked their determination of the four-loop 
RG functions by the technique of 1/n expansion. In the limit of large n, Eqs. Q2.24|) 
and ( |2.25| ) reduce to 



n" 
n 3 



2tt 2 
1 

2^ 



1 

32 
1 

32^ 



Gi 
2 



7T 



(K + 4K 2 - 2^ 



+ 



32^ ^~2^~16 + T + T + ~8~ 



+ o[- 

n 



(2.26) 
(2.27) 



Eqs. (|2.26|) and ( |2.27| ) agree with the predictions by the 1/n expansion in the given 
order JTS| . However, as we see in the above equations, this method can not control 
all integrals in Table [I], merely those of G\, K, Vi, V 2 and V 3 . 

There was, on the other hand, a further check on the RG functions when we 
determined a mass gap |T7| by applying a perturbative computation of the mass 
gap M(L) in finite volume L. In converting our computation of M(L) on the lattice 
to the renormalized MS scheme, we made there an explicit use of the four-loop /3- 
function. From the constraint that we should formally have a free theory in the 0(2) 
model with dimensional regularization, we get a relationship for the perturbative 
M(L) to satisfy. In this case, we could however confirm only the integral J. 



3 Application to physical quantities 

In this section we present the application of the perturbative determination of the 
RG functions in the previous section to the physical predictions, explicitly the an- 
alytic correction to the correlation length £ and spin susceptibility \ as well as the 
determination of a mass gap through a perturbative computation of a finite volume 
mass gap M(L). We first consider the application to the analytic correction to £ 
and x- 



3.1 Perturbative evaluation of £ and x 



The determination of the four-loop j3- and 7-functions makes it possible to compute 
£ and x U P to second analytic correction. We discuss here the numerical evaluation 
of the analytic corrections to £ and x an d also to the other relevant quantities. We 
first consider the correlation length £. 

The correlation length is a physical quantity and correspondingly satisfies the 
RG equation: 



£(/o,a) = 



df da 

We may then define the RG invariant parameter as 



(3.1) 



t/a = Q(6i/o) S2/S? exp( r -)£(/ ) 



(3.2) 



6 



£(/o) 



cxp 



Jo V /?(.-. 



+ 



o \(3(x) hx 2 b\x 



12, 



(3-3) 



which satisfies Eq. (|3.1| ) order by order in perturbation theory. The constant is 
non-perturbative and its value depends on the explicit definition of the correlation 
length. For the isovector exponential correlation length £y , which controls the 
large distance behaviour of the two-point function (q(x) ■ q{0)), an exact determina- 
tion was achieved by Hasenfratz and Niedermayer || using the thermodynamic 
Bethe ansatz: 



CV(exp) 



l/(n-2) 



r 1 + 



n 



2 -5/2 



exp 



7T 



2(n-2) 



(3.4) 



An alternative definition with respect to the isovector correlation function 
(q(x) ■ q(0)) is the second-moment correlation length 



t( 2 ) 



E x \x\ 2 (q(x)-q(0)) 



(3.5) 



There are also other definitions with respect to higher isospin channels, e.g. 
the exponential correlation length related to the isotensor correlation function 
([q(x] ■ q{0)] 2 — 1/rt) and the corresponding second-moment correlation length. 
The function £(/o) of Eq- (|3.3|) is regular at fo — and may be expanded in f Q : 



A 1 f + A 2 f 2 + O(f3), 



where 



A 1 
A 2 



b 1 k) 2 + b 2 (bl-2b 1 b 2 b 3 + blb 4 ) 
2b\ 



Inserting h,b 2 ,b 3 and 6 4 from Eqs. ( p.l8|) -( ^20|) and (|2.24|) in Eqs. 
we obtain for Ax and A 2 



(3.6) 

(3.7) 
(3.8) 
and (O), 



At 
A 2 



n 



1 2-kGi )(n - 2) H 

4tt 16 V V ; 4 



192(n - 2) 2 



7T 

24 



48 



(3n - ll) 2 + 7r(4n 2 - 19n + 17) 



(3.9) 



+6n 3 - 39n 2 + 95n - 76 



32vr 2 (n - 2) 



[7n - 16 + 27r(n - 1)] 



1 73n - 164 
96tt 2 n - 2 



C(3) 



7 



24(n - 2 
+2tt 2 GI - (3 - ir)R 
1 



[-12(n - 2){n - 3) - 12vr(3n - 5) - (3n - 11)tt 2 



2 - 2 <* - - ■ -l Ll + *(n- 3)(K - 2V 3 ) 

2An- 2 3 8 V A ; 



2V^ + (n - 2)V 2 + 4\4 - 2Vg - 2Ve - 16 W 2 



(3.10) 



Now, if we substitute for the constants in Eqs. ( pT9| ) and fl3.10| ) the values in Ta- 
ble m evaluated by us, we get for A2 an estimation which is different from that by 
Caracciolo et al.: 



A, 



1 



n-2 y 
1 

("-2) 



-0.0490 -O.OWln), (3.11) 
r (0.188 - 0.122n + 0.0405n 2 - 0.0129n 3 ). (3.12) 



It is to be remarked that if n increases, the difference with the determination by 
Caracciolo et al. becomes small, which is of course to be expected since the four-loop 
coefficient is correct in the large n limit. 

Next, we would like to compute the analytic correction to the isovector spin 
susceptibility xv = J2x{q( x ) ' ft satisfies the following RG equation: 



Xv(fo,a) = 0. 



(3.13) 



We may then define the RG invariant parameter as 

2 



Xv/a 2 = C x {bx0 2bM),b " 



x(/o) 



exp 



fa 



dx 



( -l> | — J\(/n 

Ol/O' 

2 26 2 7(x) 



Cl 



\/3(x) bxx 2 b\x (3{x) bxx 



(3.14) 
(3.15) 



Eq. ( |3.14| ) with Eq. ( |3.15| ) satisfies the RG equation ( |3.13| ) order by order in per- 
turbation theory. Contrary to Cg in Eq. ( |3.2| ), there is no exact expression for 
the non-perturbative constant C x at present. It is however known in the large n 
limit 101 Bl: 



c x = - 

x 16 



1 + -(4 + 3 7c - vr - 3 7s - 7 In 2) + o{\ 



n 



n- 



(3.16) 



where 7^ is Euler's constant and 7c — 0.4861007. 

The function x(fo) °f Eq. ( |3.15|) is well-defined at fo = and may be expanded 
in fo: 

X (fo) = 1 + Bxf + B 2 f Q + 0{fl) (3.17) 
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with 



Bo 



- 2&| + b\c 2 - bib 2 di 
b\ 

\2bl + b 1 b 2 c l -b l {2b^b 1 c 2 )] 2 
2b\ 

2b\ + SiSpi - M 2 (4fr 3 + &1C2) + b\{2b 4 - 5 3 ci + S1C3) 

2&| 



(3.18) 



(3.19) 



Inserting 61, • • • , 64 and Ci, - • • , c 3 from Eqs. (g3f)-(ggg) in Eqs. ( p7T8Q and ( ^T9|) , 
we obtain for and B 2 



#1 

B 2 



2ir(n-2) 
1 

(n - 2f 



I^n + 1) - 1 + \n 2 (" - y) - ^> - 2)Gx 

(3n - ll) 2 + — (7n 2 - 30n + 17) 
1152 v ; 96 V ' 



(3.20) 



+— (n- l)(3n 2 - 13n + 51) 
96 

1 (2ra - 1) + -^(3n 2 - 13n + 15) 



57T 



1 73n - 164 
48vr 2 ra-2 



C(3) 



6(ra - 2) 



6(3n - 4) - 6vr(4n - 5) - (3n - 11)tt 2 + 8vr 2 G 2 



-2(3 - tt)/? + --— I J - + -(n - 3)(K - 2V 3 ) 

+ n[2Vi + (n - 2)V 2 + 4V4 - 2F 5 - 2V e - 16W 2 }. 



(3.21) 



Now, if we substitute for the constants in Eqs. ( |3.20|) and (|3.21|) our values in Table [l], 
we obtain 



Bo 



n-2 
1 



(-0.1888 + 0.0626n), 



(3.22) 



(n-2) 



-(0.4191 - 0.2688n + 0.0517n 2 - 0.0108n 3 ). (3.23) 



Again, for B 2 we get a prediction which is different from that by Caracciolo et al. 

Let us also consider the ratio Xv/^ 2 which has an advantage of being possible to 
compute an additional analytic correction. Writing 



2tt 1 \ -(n-i)/(n-2) 



n-2f 



i=i 



(3.24) 



where Ry = C x /C^ is a non-perturbative universal constant, the expansion coeffi- 
cients Ci, C 2 and C 3 have the expressions 



C x 



ln-1 
An n-2 



7T-2), 



(3.25) 



9 



c 2 
c 3 



n — 1 

(n - 2) 2 
n — 1 
(n - 2) 3 



n - 2) 2 Gi (3n 2 - 23n + 31) + 

96 

[10n 3 - 79n 2 + 186n - 137) 



n — 1 2n — 3 



57T 



,(3.26) 



384 



+— — (2n - 3)(3n 2 - 26n + 37) 
ly27r 



+ r(n-l)(3n-5) 

32tt 2V a ; 



48tt 



-(3ra 3 -18n 2 + 38n-27) 



n — 1 



1 



1927T 



n - 2 



(2n 2 + 57n - 134)C(3) 



48 
+ (n- 1 
n — 7 



47T 



1 



1 



(5n - 9)vr - 2(2n - 3) 
n — 3 



16 



-(x-2y 3 : 



. __ ^ + _(_2V 4 + V 5 + V 6 + 8W 2 ) 



(3.27) 



The numerical evaluations for these coefficients with our values in Table |T] yield 
n — 1 



c 2 



Co 



n-2 
n — 1 

(n-2) 2 

n — 1 
(n - 2) 3 



-0.0908, 



;-0.0316 - 0.0120n + 0.0149n 2 ), 



(3.28) 
(3.29) 



(0.0578 - 0.0564n + 0.0273n 2 - 0.0149n 3 + 0.0041n 4 ). (3.30) 



As expected, for the expansion coefficient C3 our determination does not agree with 
the estimation by Caracciolo et al. 

Finally, we would like to consider the spin-fc susceptibility |L5f defined as 



X 



(fc) 



(3.31) 



where 



q ai ■ ■ ■ q ak - "Traces" 



(3.32) 



and the "Traces" must be such that Y" fc ai ' '"' ak is symmetric and traceless. Then, from 



the standard RG arguments one has 



X ^ = Cf > exp ( — 
x \n — 



4tt 1 



2n 1 \ -[2+fc(n+Jfc-2)]/(n-2) 



2f J \n-2f 
with a non-universal constant By including the ratio 



i=l 



(3.33) 



^(fc) / 27T 1 \ -fc(«+fc-2)/(n-2) 



t=i 



(3.34) 



10 



n 


P 








7^(4) 


3 


1.85 


89.16(21) 


0.743(2) 


0.781(2) 


0.822(2) 


3 


2.25 


1049(7) 


0.861(6) 


0.897(6) 


0.928(6) 


3 


2.60 


8569(92) 


0.901(10) 


0.934(10) 


0.958(10) 


3 


3.00 


94.6(1.6)-10 3 


0.930(16) 


0.959(16) 


0.977(17) 


4 


2.50 


34.85(9) 


0.911(2) 


0.930(2) 


0.949(2) 


4 


2.80 


86.07(37) 


0.927(4) 


0.945(4) 


0.960(4) 


8 


5.80 


33.41(8) 


0.985(2) 


0.990(2) 


0.994(2) 



Table 2: Ratio of Monte Carlo results to the theoretical perturbative predictions 
at two, three and four loops for the correlation length £ in the 0(3), 0(4) and 0(8) 
model evaluated for various values of (5. 



and its expansion coefficients 



C ( k ) = 


1 k(n 




CP = 


- 2 h \ 



-(t-2), 



+ 



k(n + k- 2) 



n-2 

one gets for the expansion coefficients 



(n-2)Gi 



1 

8^ 



96 



(3n - 14) 



(k) 



B 



(k) 



2A 1 + c[ h \ 

2A 2 + A\ + cf } + 2A 1 c[ h) 



which have the numerical values 



(k) 



n-2 
1 



0.0980 - 0.0283n + 0.0908A;(A; + n - 2)J , 
0.37866 - 0.2428n + 0.08115n 2 - 0.0257n 3 



(n - 2) 2 

+A;(A; + n- 2)(-0.0363 - 0.0187n + 0.0149n 2 ) 
+0.0041A; 2 (A; + n - 2) 2 



(3.35) 
(3.36) 



(3.37) 
(3.38) 



(3.39) 



(3.40) 



Here, our value of is not in agreement with that by Caracciolo et al. 

Now, it is interesting to compare the theory with "experiment". We would like 
to take the correlation length and compute the ratio 



(3.41) 



where £mc(/3) is the Monte Carlo value for the correlation length given as a function 
of the inverse bare coupling (3 — l// and ^(P) * s the theoretical /-loop prediction 
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from Eqs. ( |3.2j )-( |3l4"D . In Table ^] we show determinations of IZ^ for some selected 
values of n and (3 by inserting our new corrected coefficient of the four-loop (3- 
function. Compared to the old estimations by Caracciolo et al., the correction of the 
four-loop /3-function causes to reduce the four-loop contribution to the perturbative 
correlation length in all of the considered cases, where the difference is maximum 
in the 0(3) model for (3 = 1.85 and becomes small if (3 or/and n increase. It 
follows that, in every case, the four- loop contributions to the correlation length are 
smaller than Caracciolo et al. have estimated, and so the discrepancies between the 
theoretical four-loop prediction and "experiment" are larger than the estimation by 
Caracciolo et al. Explicitly, at the largest (3 values investigated, they are 2.3% for 
the 0(3) model, 4.0% for the 0(4) model and 0.6% for the 0(8) model and there 
still remains some space for higher loop contributions. 



3.2 Determination of a mass gap by finite volume method 



In our work [17|] on the determination of a mass gap m by applying the perturbative 
computation of the mass gap M(L) in finite volume L, we made an explicit use 
of the four-loop /3-function in converting our computation of M(L) regularized on 
a lattice to the renormalized MS-scheme. Since the four-loop coefficient has to be 
changed, our prediction for the mass gap has to be corrected accordingly, and so we 
would like to present here the modifications caused by the change of the /3-function. 

To the determination of the physical infinite volume mass gap m from the finite 
volume mass gap M(L), we applied two different methods: one method was the 
determination by matching the behavior of the finite volume mass gap M(L) eval- 
uated at small L with that of the mass shift So = [M(L) — m]/m known at large L 
and the other was by means of a coupling g 2 {L) running with L. We analyzed both 
methods again by inserting our new corrected value of the /^-function. 

In the case of the determination of the mass gap by applying the mass shift, we 
concluded that the method is not stable and so it is difficult to make a reasonable 
estimation for the mass gap by this method. After we substituted our new corrected 
value of the four-loop /9-function, we see substantial changes, but the whole picture 
including the final conclusion does not improve the situation. Therefore, we do 
not want to present here the evaluation of the mass gap by this method. As for 
the determination of the mass gap by means of running coupling, we however see 
interesting effects through the new corrected value and therefore we would like to 
discuss in the following the new evaluation of the mass gap by this method in detail. 

In Ref. [fT7T| , we introduced a coupling g 2 (L) running with volume L defined byQ 

g\L)=2M(L)L/(n-l), (3.42) 
where M(L) is a mass gap in finite volume which we computed to three-loop order 



2 The running coupling g 2 (L) was first defined in Ref. 
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in perturbation theory. For the coefficients of /3-function defined by 



we found 



61 



62, 

(n-l)(n-2) 
1 n — 2 r 

(n - 2) 3 xi + (n - 2) 2 X2 + (n - 2) X3 + Xi 



4 (2ti) 1 



(3.43) 

(3.44) 
(3.45) 

(3.46) 
(3.47) 



The explicit expressions for the constants \ii ' ' ' j X4 can be found in appendix B 
of Ref. [0, where the values for t 1; £ 2 and £3 in their expressions have now to be 
modified by our new values: 



n-2 



(2tt)< 



(n - 2) 2 ti + (n - 2)t 2 + *a 



wi 



ith 3 



= -1.015069687473002(1), 
t 2 = -5.9168085(1), 
t 3 = -9.0937547(1). 

By inserting these corrected numbers into the expressions for Xi, 
numerical values 



Xi 

X3 

Xi 



-1.20(1), X 2 = -3.63(1), 

23.6(1), 
-5.2123414(1). 



The A-parameter in this scheme is given by 

1 



Afv 



2 



Kg 



exp 



Jo \p( x ) hx 2 b\x 



which has the following perturbative expansion up to order 0(g 



A 



(2) 
FV 



(hg 



^-b 2 /bi e -i/(b lS 2 ) 



(3.48) 



(3.49) 
(3.50) 
(3.51) 

Xa, we get the 



(3.52) 
(3.53) 
(3.54) 



(3.55) 
(3.56) 



(3.57) 



3 In the evaluation of t\, ti and £3 we took the revised value for W2 in Ref. |35| into consideration. 
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a(3) 



A 



(4) 
FV 



V 



b\ - b\b\ + 26f 6 2 6 3 - gMgg + - 6fo 4 _ 4 
+ 2b\ 9 



(3.58) 



(3.59) 



If we insert the running couplings g 2 (L) from Table 1 of Ref. ||17|| , measured for 
given L in units of the mass gap m in 0(3) model, in Eqs. ( p.57| )-( p39|) , we obtain 
new values for the ratio Apy/m in 2-, 3- and 4-loop approximations which are listed 
in Table ^] and also illustrated in Fig. [I]. 

The new effect which mainly arises by substituting our corrected value for the 
four-loop /3-function is that, compared to the old diagram with the value by Carac- 
ciolo et al., for large couplings the new four-loop estimations for Kpyjm become 
smaller, so that the four- loop curve in Fig. p] shifts underneath. As a result, the 
order of 3- and 4-loop curves changes. However, when the coupling goes to zero, the 
difference with the old diagram becomes smaller and smaller, and at the smallest 
coupling considered, it is practically zero. From this observation we can conclude 
that the changes in /3-function influence only at large couplings and for small cou- 
plings, which are indeed the region of our main interest, we have the same picture 
as the earlier investigation with the value by Caracciolo et al. Therefore, as in the 
case with the old four-loop /3-function, our determination of the mass gap converges 
to the value predicted by Hasenfratz and Niedermayer very well. 

To give a further impression on the systematic errors, we also consider an al- 
ternative definition of n-loop approximations to Apv] instead of expanding X(g 2 ) of 
Eq. (|3.56|) in the coupling g 2 , we insert the perturbative coefficients of f3(x) in \{g 2 ) 
and integrate this exactly. If we apply the coefficients of the /5-function calculated 
up to 4-loop order, we obtain the following approximations to Kpy at 2-, 3- and 
4-loop levels which are a little modified from Eqs. (|3.57|) - (|3.59|) : 



A 



(2) 
FV 



A 



(3) 
FV 



A 



(4) 
FV 



1 



L 
■ exp 



dx 



I 1 b 2 



o \x 2 (bi + b 2 x) bix 2 b\x 



1 

T 



-^i# 2 r fe2/fe? e- i/(6i52) 



• exp 



dx 



1 

T 



±r% 2 )- fe 2/fc? e -v(M 2 ) 



■ exp 



dx 



o \x 2 (bi + b 2 x + b 3 x 2 ) hx 2 b\ 



x 



+ 



o \x 2 (b 1 + b 2 x + b 3 x 2 + b 4 x 3 ) hx 2 b\ 



x , 



(3.60) 



(3.61) 



. (3.62) 
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f{L) 




A%/m 


A [ Fv/m 


0.5372 


0.0511(8) 


0.0467(6) 


0.0463(17) 


0.5747 


0.0512(7) 


0.0465(6) 


0.0460(17) 


0.6060 


0.0514(7) 


0.0464(6) 


0.0459(15) 


0.6553 


0.0518(6) 


0.0464(5) 


0.0457(15) 


0.6970 


0.0520(6) 


0.0462(5) 


0.0455(15) 


0.7383 


0.0525(6) 


0.0463(5) 


0.0455(13) 


0.7646 


0.0526(6) 


0.0462(5) 


0.0453(13) 


0.8166 


0.0532(6) 


0.0463(5) 


0.0453(11) 


0.9176 


0.0552(5) 


0.0471(5) 


0.0458(10) 


1.0595 


0.0574(4) 


0.0478(3) 


0.0459(8) 


1.2680 


0.0628(1) 


0.0502(2) 


0.0473(8) 



Table 3: Two-, three- and four-loop approximations to Kpyjm m the 0(3) model 



0.08 



0.07 



B 

\ 0.06 



0.05 



0.04 

1.2 1 0.8 0.6 

g a (L) 



Figure 1: Two-, three- and four- loop approximations to A FV /m in the 0(3) model 
from Table |3|. The triangular points show the 2-loop approximation to the A- 
parameter, the hexagonal points the 3-loop and the circular points the 4-loop approx- 
imation. The dotted line is the value of Hasenfratz et al. {Kpy /m = 0.04815 ■ ■ ■). 
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Insertion of the data in Table 1 of Ref. |I7| in these equations gives the numbers 
in Table |] which are also plotted in Fig. |2|. In this case, we have exactly the same 
phenomenon as the former investigation. Compared to the old figure with the four- 
loop /3-function value by Caracciolo et al., the four-loop curve for Apy/m at large 
couplings shifts underneath in this well, so that the splitting of 3- and 4- 

loop curves becomes larger. However, for small couplings which are actually our 
interesting region, the difference between two figures, one figure with the new (3- 
function value by us and the other with the old value by Caracciolo et al., vanishes. 
Therefore, also in this investigation with the modified expansion to the A-parameter, 
we can conclude in the same way as the earlier investigation with the old value of 
the /3-function: all three curves converge to the value predicted by Hasenfratz et al. 

We also want to discuss the changes of the diagrams through the correction of 
the /^-function in 0(4) model. By substituting the new value of the /3-function in 
Eqs. Q3.57) )-( pr5S| ), we obtain the corrected results for A FV /m which are written in 
Table [| and also illustrated in Fig. [| As far as the modification of the four-loop 
curve is concerned, Fig. ||] shows exactly the same behavior as Fig. [I]: as a result of the 
shift of the four-loop curve underneath, the order of 3- and 4-loop curves changes, 
but for small couplings the difference between the old diagram and the new one 
vanishes and we have the same argument with respect to the final conclusion. 

As for the determination of Apy/m by means of the modified definition of the 
perturbative Apy in 0(4) model, we obtain the results which are written in Table |6] 
and also plotted in Fig. f|. In this case, we see exactly the same phenomenon as 
Fig. ^| which leads to the conclusion that for small couplings all three curves converge 
to the value by Hasenfratz et al. very well. 

Summarizing the whole discussions up to now, we conclude that the correction 
of the four-loop /5-function has an influence on our earlier investigation for the de- 
termination of the mass gap only at large couplings, while in the actually interesting 
region of small couplings it provides practically no changes. Therefore, our conclu- 
sions with the old value of the /3-function by Caracciolo et al. remain valid and we 
can say that in all of the investigated cases our determinations of the mass gap with 
the corrected value of the four-loop /3-function by us agree with those by Hasenfratz 
et al. very well. 
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9 2 (L) 


Afy/m 


Kf v /m 




0.5372 


0.0555(8) 


0.0474(6) 


0.0465(17) 


0.5747 


0.0559(7) 


0.0473(6) 


0.0462(17) 


0.6060 


0.0563(7) 


0.0473(6) 


0.0461(15) 


0.6553 


0.0572(6) 


0.0474(5) 


0.0460(15) 


0.6970 


0.0578(6) 


0.0474(5) 


0.0459(15) 


0.7383 


0.0586(6) 


0.0476(5) 


0.0459(23) 


0.7646 


0.0590(6) 


0.0476(5) 


0.0458(13) 


0.8166 


0.0601(6) 


0.0479(5) 


0.0458(11) 


0.9176 


0.0633(5) 


0.0492(5) 


0.0466(10) 


1.0595 


0.0671(4) 


0.0505(3) 


0.0472(8) 


1.2680 


0.0755(1) 


0.0544(2) 


0.0496(8) 



Table 4: Alternative 2-, 3- and 4-loop approximations to Kpyjm in the 0(3) model 



0.08 n 



0.07 



"i 0.06 



0.05 



5 i J - 



GO O O O O 



0.04 



1 .2 
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0.6 



I 2 (L) 



Figure 2: Two-, three- and four- loop approximations to A^y/m in the 0(3) model 
from Table 0. The arrangement of the points and the dotted line is the same as in 
Fig. 0. 
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Kfy/m 


Kf v /m 




0.863 


0.0795(7) 


0.0722(7) 


0.0717(16) 


1.011 


0.0817(7) 


0.0729(6) 


0.0722(15) 


1.228 


0.0844(5) 


0.0734(6) 


0.0723(12) 


1.584 


0.0880(5) 


0.0733(6) 


0.0713(10) 


2.309 


0.0928(4) 


0.0701(4) 


0.0657(9) 


4.132 


0.0853(4) 


0.0479(4) 


0.0351(9) 



Table 5: Two-, three- and four- loop approximations to Apv/ m m the 0(4) model 



0.1 



0.08 - 



< 



0.06 



0.04 




Figure 3: Two-, three- and four- loop approximations to Apy/m in the 0(4) model 
from Table |5|. The arrangement of the points is the same as in Fig. |l|. The dotted 
line is the value of Hasenfratz et al. (Apv / m — 0.07321 ■ ■ •). 
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0.863 


0.0830(7) 


0.0733(7) 


0.0720(16) 


1.011 


0.0859(7) 


0.0744(6) 


0.0726(15) 


1.228 


0.0898(5) 


0.0757(6) 


0.0731(12) 


1.584 


0.0952(5) 


0.0771(6) 


0.0730(10) 


2.309 


0.1035(4) 


0.0780(4) 


0.0709(9) 


4.132 


0.1024(4) 


0.0678(4) 


0.0561(9) 



Table 6: Alternative 2-, 3- and 4-loop approximations to hpylm in the 0(4) model 



0.1 



0.08 



0.06 - 




I 2 (L) 



Figure 4: Two-, three- and four-loop approximations to Kpyjm in the 0(4) model 
from Table || The arrangement of the points and the dotted line is the same as in 
Fig. H 
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A Evaluation of lattice Feynman diagrams 

In this appendix we define the lattice integrals appearing in Table [1] of section ^| and 
show how they can be evaluated. 

A.l Evaluation of W 2 

W2 is defined as a finite part of an infrared divergent integral: 



W 2 = lim 
with 



TF 2 (,)->)3 + (i-^)^) 2 -(l + ^- T i--f)^) 



(A.l) 



jt , M / E^sin^sinrvsin^sin^ 1 - 

W ^ ) = I (? + V<P + ^ + V<P + V {q 7rf + h Siq + r + 3 + th <A - 2) 



qrst 

\2; 



where 5(g) stands for a modified form of 5-function (2n) 5(q) and / for an abbrevi- 
ated symbol of the integral sign J^ n t^w- The quantity R is a finite integral which 
can be evaluated exactly 



R = lim 



1 (A - 3) 



24vr 2 



with = dlnr(z)/d2;, while the function if(/i) is a divergent integral for small h 

w) = fpTH = -h ,n m +0 <- hblh) (A ' 4) 

p 

with p 2 = J^fiPfi, fin = 2 sin(p At /2) and /i = 0, 1. Our aim is to evaluate the integral 
W2{h) in the limit h — > to confirm the logarithmic terms in Eq. (|A.1|) and to 
determine the constant W2 a t the end. 

Observing that Eq. (|A.2|) can be represented by the following symmetric form 

WW= / ^rE^fc)' (A.5) 

with 

w^Z ^^V ^^ (A - 6) 

we try to evaluate Eq. (|A.5|) by dividing into the following two integrals: 

W 2 (h)=h(h)+I 2 (h), (A.7) 
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where 



h(h) 



P 2 + h ^u 



(A.8) 
(A.9) 



fJ,,U 



Considering, at first, the integral h(h), we transform the function 0) to 

position space to study how it behaves for small p: 



where H(x) is defined by 



H(x) = J e ipx \np 2 . 



(A.10) 



(All) 



For large x, Eq. ( JA.ll ) shows the behavior lim H(x) ~ js, from which we can find 
out that F^lp, 0) diverges logarithmically for p — > 0. We then break up Eq. ( |A.1(J| ) 
into converging and diverging part 



o) + J- £ e^(9, + + flj)c?(x) - E e^to 

47T ~ r Z7T ~ 



(A.12) 

where the free propagator including its properties is discussed in detail in 

Ref. HH and F^p, 0) is a well-defined function in the limit p — > 0, constructed by 
subtracting from Eq. ( |A.1U| ) the leading square terms in large x: 



f»Ap, o) = £e^ ^h(x) 2 - —(d, + o;)(d u + d;)G(x) + —s^h 



2tt 



x , 



(A.13) 



Next, we separate the p = part of Eq. ( |A.13j ) to introduce a function F^ip, 0) 
which vanishes for p = 0: 



F„ u (p, 0) = F„ u (p,0) + 6 fMU C 



(A.14) 



with 



F^{p, 0) = - 1) U^H(x) 2 - -d,d;G(x) + S ^H(x) 

l 71 Lix 



, (A.15) 
(A.16) 

The numerical value of C, which can be computed very precisely by applying the 
position space method |L6|], is 0.369899 • • •. 
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Substituting Eq. ( |A.12| ), together with Eq. ( |A.14| ), in Eq. ( |A.8| ) and going back 
to momentum space, we have 



h(h) = - 



l r l 



16 7 p 2 + h 



E F^M 



1 sinp^sinp,, 1 ^ 
b — V ln P + <W 

71 p z Z7T 



(A.17) 



with 



*MP,0) = -4 / ^- S(p+q+r) + - 

J q z r z 7i 



1 sin p M ship,, 1 . 2 

2 h— Vlnp ( A -^ 



2tt 



Expanding the square, Eq. (|Q7D prod uces the following ten integrals which we can 
evaluate separately: 



2C 2 A'(A) - -C/f(A) + ^/f(A) - -C/f(fc) + -U J (ft) 



71 



2n 2 ' 



7T 



7T 



1 



-Jt{h) + 2C//(/i) - -Jf(/i) - -It(h) + 11(h) 



7T Z 



71 



TV 



(A.19) 



with 



11(h) 

11(h) 
If(h) 

11(h) 
I{(h) 
11(h) 
iUh) 
I[(h) 



Inp 2 

p 2 + h' 



(Inp 



2\2 



p 2 + h 



1 E M sin p M 
p 2 + h p 2 

lnp 2 E^sin 2 ^^ 
p 2 + h p 2 

1 ^ sin 2 p M sin 2 p v 

2\2 ' 



E 



P 2 + hrr^ (p 2 ) 



1 



p 2 + h ^ 

Inp 2 

p 2 + h (l 



E^feO), 

A* 

E^>>°)> 



p 2 + h 



smp^ sinpj, 



(A.20) 
(A.21) 
(A.22) 
(A.23) 
(A.24) 
(A.25) 
(A.26) 
(A.27) 
(A.28) 
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The first five integrals diverge for small h. It is not difficult to first compute the 
divergent terms analytically and then evaluate the remaining constant terms numer- 
ically in the limit h — > 0. Our results are 



11(h) 

11(h) 

I((h) 
I{(h) 



-^-(\nh) 2 + Cl + 0(h\nh), 

—\(lnh) 3 + n 2 lnh\ + c 2 + 0(h\nh), 

K(h)-\c 3 + 0(h\nh), 

7T Z 

-— (Inh) 2 + c 4 + O(hlnh), 

8n 

K(h) + c 5 + 0(h\nh), 



(A.29) 

(A.30) 

(A.31) 

(A.32) 
(A.33) 



where the values for c\ , • ■ ■ , C5 are listed in Table [7]. The last four integrals in 
Eqs. ( |05D -( [7C28D are, on the other hand, well-defined for h -> 0. Eqs. (|AT25| )- 
( A. 27 ) can be transformed to lattice sums in this limit 



7/(0) = Y.G(x)\x 2 H(xf-^-Y.^ + d;fG(x) + lH(x) 

7f(0) = Y, L (x)\x 2 H(x) 2 -^Y,(d, + d;) 2 G(x) + -H(x) 

C 
+ 4' 



I*(0) 



~ ^G(x)^x I/ [y Ali/ (x + /i) - V^(z-/z) 



(A.34) 
(A.35) 
(A.36) 



where the functions and V^ u (x) are defined by 

hip 2 



Ux) 



(A.37) 



^(x) = x,x u H(x) 2 -^(d, + d;)(d u + dt)G(x) + ^5, u H(x). (A.38) 



The numerical values for Eqs. ( |A.34|) -( [A.36| ) which we computed by applying the 
coordinate space technique are given in Table [7]. As for the last integral I[(0) of 
Eq. (|A.28 ), we can compute it by a Fortran program calling NAGLIB integration 
subroutines. We get a number in Table |7|. 

Let us now consider the integral I 2 (h) of Eq. ( |A.9| ). With help of the function 



H^ip, h) = F^ u (p, h) - F^(p, 0), 
it can be divided into two integrals 



(A.39) 



I 2 (h)=A(h) + 2B(h) 



(A.40) 
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with 



A(h) 
B{h) 



p 2 + h 



p 2 + h 



(A.41) 
(A.42) 



which we would like to evaluate in the h — » limit. 

Noting that Eq. ( |A.41 ) has a finite limit for h — > 0, we perform the variable 
change p — > Vhq and take h — > to transform to 



A(Zi) 



h^o r d g :i 

loo (27T)V + lfe ""^ 



(A.43) 



with 



= hm H^(qVh,h) 



h-+0 



oo cl 2 ^ 



3 (27r) 2?>S "Vr 2 + ls 2 + l r 2 s 2 
which can exactly be evaluated by applying Feynman parametrization 



q — r 



(A.44) 



-5^1 -lng 2 + 



2 V? 2 + 4 arctanh-i / 



<T+4 



arctanh v /^ 



71 



(q^sfq^TA 



. (A.45) 



By transforming to polar coordinates, the two-dimensional integral A(0) of 
Eq. (|A.4ij[ ) can be reduced to an one- dimensional one being available to a numerical 
integration to high precision. Its result is reported in Table [T[ 

The integral B(h) from Eq. QA.42 ) is a little more tricky since it diverges for 
h — > 0. We apply Eq. ( |A.12j ) to rewrite 



B(h) 



1 r 



B 1 {h)+B 2 {h) 



where 



B 2 {h) 



«-rE^(P.^{P.°). 



f + h Jj 



J^H^ip, h) 



Inp 2 — — 

2lT 7T 



1 sin p,, sin p,, 
~ r 



(A.46) 

(A.47) 
.(A.48) 



For the evaluation of Bi(h), we perform the variable change p — > and let h — > 0. 
Then, we have 
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Constants 


Values 


Constants 


Values 


Cl 


0.150666215804(1) 


if(0) 


-0.0165452709(1) 


c 2 


0.420051822207(1) 


iftO) 


0.0169938575(1) 


c 3 


1.682002936874882(1) 


A(0) 


0.001687475879587(1) 


c 4 


-0.0355227164897(1) 


B« 


0.017838149659537(1) 


c 5 


-0.2327464349513382(1) 




-0.006827009250335(1) 




-0.066086583814(1) 




0.001497591836898(1) 


/f(0) 


-0.07360068(1) 


5 2 d 


0.017838149659537(1) 



Table 7: Numerical values for constant terms appearing in the determination of 
W 2 . 



which is identically zero. As for the integral B2(h), we do the same variable change 
p — > \/hq. In the limit h — > 0, we get 

B 2 (h) h =° -±-{\nh)B a 2 - ±B b 2 - -B c 2 + CB d 2 (A.50) 

Z7T ZTX 71 

with the two-dimensional, well-defined integrals 

£<§vti (a - 51) 

d 2 q Inq 2 



B« 



B\ 



B c 2 



B d 2 



Jl^fkV M (A - 52) 



]2 



/I^^E^W^. (A. 53) 

,$r£lX B M- (A.54) 



By transforming to polar coordinates, all of these integrals can again be reduced to 
one-dimensional ones whose numerical values are listed in Table [7|. 



A. 2 Evaluation of W\ 

W\ is defined as a finite part of an infrared divergent integral: 



W\ = lim 



Wl (h) + ^K(hf-^K(h) 



with 



W 1 {h) 



qrs 



(q 2 + h) 2 (f 2 + h)(s 2 + h) 



S(q + r + s). 



(A.55) 



(A.56) 
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Our aim is again to evaluate the integral W\{h) in the limit h — > to confirm the 
logarithmic terms in Eq. ( A.55Q and to determine the constant W\ at the end. 
Eq. ( |A.56p can be transformed to a symmetric form 



qrs 



(q 2 + h) 2 (f 2 + h)(s 2 + h) 



5(q + r + s) 



(A.57) 



which can, by expanding the numerator, be divided into two integrals 



W 1 (h) = 2T 1 (h)--T 2 (h) 



with 



Ti(h) 



qrs 



qrs 



fi sin s M 

(? 2 + /i) 2 (f 2 + /i)(s 2 + /i) 

{q 2 + h) 2 {f 2 + h){s 2 + h) 



5(q + r + s) 



£(g + r + s). 



(A.58) 

(A.59) 
(A.60) 



Noting that Eq. (|A.59|) can, with help of Eq. (|A.6| ), be represented by 



(A.61) 



(q 2 + h) 2 ' 

H 

we can evaluate it analogously to Eq. ( |A.5| ). We break up Eq. ( |A.61| ) into two parts 



T 1 (h)=T?(h) + T b l (h), 



where 



T*(h) 



E 
E 



(q 2 + h) 



2 F w(<li °)' 



F^(q, h)-F^(q,0) 



{q 2 + h) 2 

By analogy to h{h) of Eq. ( |A.8| ), Eq. ( |A.63| ) can be expanded to 



(A.62) 

(A.63) 
(A.64) 



W) = -iE 



4V (Q 2 + h) 2 



1 sin g M 1 



7r g z 



— (J^ ln<f + C 
In 



(A.65) 



which produces the following four integrals being able to be evaluated separately: 



Tf{h) 



CJi(h) - ^-J 2 (h) - -J 3 (h) + J 4 (h) 



(A.66) 
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with 



Ji{h) 



(q 2 + h) 2 ' 

,2 



Q 



(q 2 + hj- 



lnq 2 



E 
E 



ql sin 2 
(q 2 + h) 2 q 2 ' 



2 ^(?,0). 



(q 2 + h) 



(A.67) 
(A.68) 
(A.69) 
(A.70) 



The first three integrals diverge for h — *■ 0. After having isolated the divergent parts 
analytically, the remaining constant terms can be determined numerically: 



Uh) 
Uh) 

Uh) 



K(h) - — + 0(h\nh), 

4:71 



16n 



(\nh) 2 + 2\nh +k! + 0(hkih), 
\nh + k 2 + 0{h\nh), 



(A.71) 
(A.72) 
(A.73) 



where the values for k\ and hi are given in Table [8]. The last integral J^{h) is, on the 
other hand, infrared finite and can, in the h — > limit, be transformed to a lattice 



sum 



U$) = -EE [(d,-d;)G 2 (x) 



X fl 



xlH{xf - -L(0„ + d;fG{x) + ±-H{x) 



(A.74) 



where the function G 2 (i) including its properties is discussed in detail in Ref. |16 
The number of J^O) computed by the position space technique is shown in Table | 
Concerning the integral T^(h) of Eq. ( |A.64|) , it can be expressed by 



(A.75) 



with Hpufa, h) from Eq. (|A.39|) . By changing the variable q = \fhp and taking the 
limit h — > 0, Eq. ( |A.75| ) becomes 



2*(0) = £ 



oo ({ 2 p 



(2tt) 2 (p 2 + 1) 



2 h Up)- 



(A.76) 



By transforming to polar coordinates, this two-dimensional integral can be reduced 
to one-dimensional one whose numerical value is reported in Table M. 
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Constants 


Values 


Constants 


Values 


ki 


0.150666215804(1) 


J 4 (o) 


-0.057000592735(1) 


k 2 


0.00658179485930(1) 


I*(0) 


0.001497591836898(1) 


h 


0.68169011381621(1) 


IJ(0) 


-0.02533029591(1) 



Table 8: Numerical values for constant terms appearing in the determination of 



The last integral T 2 (h) of Eq. ( |A.60| ) is to be evaluated by dividing into two parts 

T 2 (h)=T 2 a (h)+T%(h), (A.77) 

where 



T?(h) 



qrs 



(q 2 + h) 2 (f 2 + h)(s 2 + h) 



5(q + r + s)- S(r + s) , (A. 78) 



(f 2 + h)(s 2 + h) 



5(r + s) 



(A.79) 



In the limit h — > 0, Eq. ( |A.78|) is finite and can be transformed to a lattice sum 



T 2 a (0) = - EE [(^ - 9;)G 2 (x)} [(0„ - d;)G(x) 



(A.80) 



X (J, 



while Eq. ( |A.79| ) reduces to 



T b 2 (h) = 



+ 0(Mn/i) 



with 



y (^4. m 



(A.81) 
(A.82) 



The numerical values for T 2 a (0) and k% are listed in Table 



A. 3 Evaluation of the constant lattice integrals 

The constant lattice integrals appearing in Table |l]of section [| are defined as follows: 



J 



1 M 2_r 2 -s 2 E»qt 



f 2 s 2 



qrs 



2^2 



8(q + r + s), 



qrst 



q 2 f 2 s 2 t 2 



-S(q + r + s +t), 
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K 

Li 
Vi 

V-2 



- — - 2 i r - — - 2 

(q + r) — q 2 — f 2 (s + 1) — s 2 — t 2 



qrst 



q 2 r 2 s t 2 



S(q + r + s + t), 



A. S\ A " " r 

T,fM,v q^qvS^Syt^ty r^r v (s + t) - (s + t)^(s + t)„r* 



qrst 



qrst 



^2fi q^ r fi s fi^/j, 

q 2 r 2 s 2 t 2 



(q 2 ) 2 r 2 s 2 i 2 (s + tY 
5(q + r + s + t), 



5{q + r + s + t), 



s + tf - s 2 - t 2 ] [\(q + r) 2 <f + (qTr)X 



qrst 



S 2 t 2 



E 



(q 2 ) 2 r 2 



x 



+ rl)(q + r) 2 - (q 2 + f 2 )(q + rf 



[(q + r) 



212 



- — - 2 ~ i r - — - 2 

(s + t) - s 2 - i 2 (q + r) - q 2 - f 2 



qrst 



q 2 f 2 s 2 t 2 [(q + r) 2 ]' 



5(q + r + s + t), 



qrst 



q 2 f 2 s 2 t 2 (q + r) 2 



S(q + r + s + t), 



qrst 



J2^,u q^s^ sin 2 (g + r) u - 
q 2 f 2 s 2 t 2 (q + r) 2 



S(q + r + s + t), 



[s + t) - s 2 - i 2 ] J2fi (q + r) 2 ^^ cos \(q + 



IP r 



qrst 



q 2 f 2 s 2 i 2 (q + r) 



S(q + r + s + t). 



We can transform all of these integrals except V4 and V5 to position space and 
then apply the coordinate space method ||16|| : 

G i = lE{T,d;d,d;d,G 2 (x)}\Y,[((d,~d;)G(x)) 2 + ((d, + d;)G(x) 

X f£ ' fjb 

X fl,V 



K = \E{J2[n-9;)G(x)) 2 + ((d ll + d;)G(x: 



d^d u G(x) 



Li = EE[W 

X fl,U 

vi = EE[W 

X fl 



i 2 



x v /J, 



+ J2[(d, + d;)G( 



x) 
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X 



G(x) - d;fG 2 {x) + £ [(d, - d;)G(x)} [(d, - d;)G 2 (x) 

/X fj, 

v, = Eg(*) 2 E(^-^) 2 g(s) 2 , 

X fl 

V 6 = WG(x) 2 J2{2[(d,-d;)(d, + d;)G(x)}[G(x + 2^-G(x-2^ 

" x n 

+ [{d, + d;fG{x)] [{d, - d;fG{x) + (d, + d;fG{x)] }. 

Concerning the integral V4, after transforming to a symmetric form 

V 4 = -J iEVW 2 (A.83) 

with 



f 2 

t 1 w 



(t + q)l sin ^ 



(A.84) 



we can compute it by a Fortran program calling NAGLIB integration subroutines. 
The integral V5 can, on the other hand, be divided into 



where 



v 5 a = 



qrst 



. 5 ( q + r + s + t ), 
q A r s t 



^ / 



q 2 f 2 s 2 i 2 (q + r) 2 
V 5 a is then able to be represented in coordinate space 



(A.85) 

(A.86) 
(A.87) 



(A.88) 



x 11 



while V 5 6 can be symmetrized as 



v 5 b = j^J2BM 



with 



B v (t) 



(t + q) v q v 



(t + qW 

which is again to be evaluated by using NAGLIB integration subroutines. 



(A.89) 



(A.90) 
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